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ABSTRACT 

Primordial stars formed in the early universe are thought to be hosted by compact dark matter 
(DM) halos. If DM consists of Weakly Interacting Massive Particles (WIMPs), such stars may 
be powered by DM annihilation during the early phases of their evolutions. We study the pre- 
main sequence evolutions of the primordial star using a detailed stellar evolution code under the 
assumption that the annihilation of adiabatically contracted WIMPs DM within the star provides 
a sufficient energy to sustain the stellar equilibrium. We follow the evolution of accreting stars 
using several gas mass accretion rates derived from cosmological simulations. We show that the 
stellar mass becomes very large, up to 900 — 1000 Mq when the star reaches the main-sequence 
phase for a reasonable set of model parameters such as DM particle mass and the annihilation 
cross section. During the dark star phase, the star expands over a thousand solar-radii, while 
the surface temperature remains below 10^ K. The energy generated by nuclear reactions is not 
dominant during this phase. We also study models with different gas mass accretion rates and 
the DM particle masses. All our models for different DM particle masses pass the dark star 
phase. The final mass of the dark stars is essentially unchanged for DM mass of m^^ < 10 GeV. 
Gravitational collapse of the massive dark stars will leave massive black holes with mass as large 
as 1000 Mq in the early universe. 

Subject headings: cosmology: first stars — cosmology: dark matter — stars: evolution 

1. Introduction ture observations of the distant universe will ex- 

ploit large ground-based and space-borne tele- 
The first stars in the universe may have con- g^^^^^ ^^^^i as James Webb Space Telescope and 

tributed early cosmic reionization and may have j^^i^^y y^^^^j. Telescope to answer the important 

also enriched the inter-galactic medium with questions of how and when the first stars were 

heavy elements such as carbon, oxygen and iron formed and how they affected the subsequent evo- 



see, e.g. iBromm et al.l (|2009l ) for a review). Fu- lution of the universe 
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Recent theoretical studies based on cosmolog- 
ical simulations suggest that the first stars are 
formed in dark matter (DM) h a.los with mass 
10^ - 10^ Mr at z - 20 - 30 (IXegmark et al 



1997t lYoshida et al.ll2003h . In such a cosmological 
"minihalo" , the gas cools and condenses by molec- 
ular hydrogen cooling to form a star-forming gas 
cloud. The gravitationally unstable cloud further 
co ntracts, and finally a proto-star is born in side 



it ()Omukai fc Nishilll998t lYoshida et alJl2008l) . A 
unique characteristic of primordial star formation 
in the standard cold dark matter model is that the 
star and its parent gas cloud are embedded at the 
center of the host dark halo. Thus the formation 
and evolution of primordial stars may be much af- 
fected by DM dynamically. 

The nature of DM remains still unknown, but 
the best candidates are thought to be Weakly 
Interacting Massive Particles, WIMPs. WIMP 
DM such as neutralinos must have a large self- 
annihilation cross-section in order to be a dom- 
inant component of DM in the present-day uni- 
verse. DM annihilation produces an enormous 
amount of energy, essentially equal to the rest 
mass energy of the annihilated DM particles. If 
DM density can be very large in a primordial gas 
cloud and in a primordial star, the annihilation 
energy would act as an efficient heat input. DM 
annihilation could even supply sufficient energy to 
support self-gravitation of a star. This is indeed 
a recently proposed new stellar phase, in which 
stars are powered by DM annihilation ener gy in- 



stead of that of the nuclear fusion ( Spolvar et al 



20081 ). Hereafter we call such a star "dark star" 



Spolvar et al.l (|2008[ ) suggest that, in a collaps- 
ing primordial gas cloud, DM annihilation heating 
can win over the radiative cooling , effectively halts 



further collapse of the cloud. iNataraian et al 
(|2009[ ) study the evolution of the DM density pro- 
file during the first star formation using cosmolog- 
ical simulations. They show that the DM density 
profile indeed becomes very steep, as much as ex- 
pected by adiabatic contraction models. 

There have been many previous studies on dark 
stars. While most of them assume constant stel- 
lar mass models ( locco et al.l 2008 : Taoso et al 



2008 


; Yoon et al. 


and 


Umeda et al. 




Spolvar et al.l (l2009^ 



study dark star evolu- 
tion with gas mass accretion. By using an ap- 



proximated stellar structure model, ISpoIv ar et al 



()2009l) follow dark star evolutions with adiabati- 
cally contracted DM (ACDM). They show that the 
final stellar mass be comes as large as 1000 Mq. 
lUmeda et al. (l2009l) perform detailed, dynami- 
cally self-consistent stellar evolution calculations 
up to gravitational core-collapse. They show that 
the final mass can be as large as 10^~^ Mq or 
more, if the star captures DM efficiently. However, 
the DM capture rate itself is uncertain because it 
depends on the unknown DM-baryon scattering 
cross section and also on the ambient DM density 
during the dark star evolution. Thus the proposed 
formation of inter-mediate massive black holes ap- 
pears to occur only in very particular cases. 

An attractive feature of the dark star models is 
that the possibility for the birth of massive stars, 
and black holes (BHs) as the end of such stars, in 
the early universe. Super-massive BHs (SMBHs) 
existed already at z ^ 6, when the age of the uni- 
verse is less than one billion years old, but the for- 
mation mechanism of such early SMBHs remains 
unknown. Gravitational collapse of a very massive 
dark star might provide a solution to this prob- 
lem. We explore a possibility that the first stars 
formed in DM halos can become very massive by 
being powered by DM annihilation. To this end, 
it is important to determine the final state of the 
dark star phase exactly, and to determine the final 
stellar mass. 

In this paper, we study the evolution of an ac- 
creting dark st ar using a detailed stellar evolu- 
tion code as in Umeda et al. 1 (l2009l) but with the 
ACDM annihilation. We calculate the density pro- 
file of ACDM using an analytical method and de- 
termine the DM annihilation rate. We then follow 
the evolution of the dark star which grows in mass 
by gas accretion. The dark star model has two pa- 
rameters effectively; the gas mass accretion rate 
and DM particle mass. We calculate several mod- 
els with different set of parameters and investigate 
the effects on the final stellar mass. Especially, we 
clarify the effect of the different mass accretion 
rates, w hich were not explicit ly shown in a similar 
work by ISpolvar et cH (|2009[ ) . 

The rest of the paper is organized as follows. In 
section 2, we introduce our numerical calculation 
and explain how we implement DM annihilation 
in the stellar evolution. Sec. 3 shows the results 
of our dark star models and we discuss the impli- 
cations in Sec. 4. We give our concluding remarks 
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in Sec. 5. 

2. Method 

We use a stellar evolu tion code with gas 
mass a ccretion devel o ped b y lOhkubo et aL I (I2OO6I 



2009t) : lUmeda et all pOOOl ). We implement en 



ergy generation by DM annihilation rate in the 
code in the following manner. We consider a 
primordial (proto-)star embedded at the cen- 
ter of a small mass DM halos. The DM den- 
sity profile, Pv, is calculated usi n g the analyt- 



ical method of Blumenthal et al. (|l986 h: when 



the primordial gas cools and collapses, DM par- 
ticles are also "dragged" into the gravitational 
potential well. The adiabatic contraction model 
uses an adiabatic invariant. For spherically dis- 
tributed DM and baryons with a total mass of 
M{R) — Mgas(-R) + AI-dm{R), an adiabatic invari- 
ant is 

M{R)R^ const. (1) 

Freese et al. I (l2009l) compared the difference of DM 
density profile calculated using this method and a 
more accurate method developed bv lYoung (1980). 
They conclude that the difference is not more than 
a factor of two. 

We have performed direct numerical simula- 
tions of early structure formation to check the 
accuracy of this analytical method. To this end 
we have used the parallel iV-body/Smoothed Par- 
ticle Hydrodyn amics (SPH) solver GADGET-2 
()Springell l2005l) m its version suitably adopted 
to follow radiative cooling pro cesses at very high 



densities (jYoshida et al.l 120061 ). We have found 



that the DM density calculated by Blumenthal's 
method is indeed in good agreement with the re- 
sult of our direct numerical simulations. Figured] 
compares the results from the analytic model and 
the simulations. Details will be presented else- 
where (Hirano et al., in preparation). We note 
that, although these two results agree very well, 
confirming the validity of the analytic model, the 
simulations do not resolve the DM density down 
to length scales of astronomical units. Because the 
gas dens ity profile itself evolve s in a self-similar 
manner (jOmukai fc Nishil 119981 ) . we assume that 
the DM continue to contract adiabatically and the 
evolution of DM density can be calculated by Blu- 
menthal's method. 

In this work, we only consider the DM den- 



sity evolution based on adiabatic contraction. In 
principle, DM annihilation can occur even more 
efficiently in the star by capture. If the den- 
sity of the stellar core increases enough to scat- 
ter and trap DM particles, the DM particles can 
rapidly sink toward the center, to self-annihilate 
efficiently. During the dark star phase, however, 
the effect of DM capture is unimportant because 
the star has expanded extremely and the gas den- 
sity is low. 

Suppose a "cloud" of DM particles contract 
from a radial position i?pro to i?now, and the to- 
tal mass distribution changes from M-pi-c{R) to 
M^cvj{R)- Then the following relation holds dur- 
ing adiabatic contraction: 



(2) 



The "new" DM density profile can be simply cal- 
culated from the above equation. Note that, in our 
case, the gas distribution also changes because of 
the gas accretion onto the central primordial pro- 
tostar. We model the gas mass accretion rate in 
a simple parameterized form as a function of the 
stellar mass. We use the accretion rate calculated 



by cosmological simulations of iGao et al.l (|2007l ). 
They found a large variation of accretion rates. 
In particular, the accretion rates for rotationally 
supported disks are typically smaller than for the 
other cases. As our fiducial model, we choose the 
accretion rate of their R5 run, which is well de- 
scribed as a power-law 

dM 



dT 



0.18 X A/-°-^ Mo yrs" 



(3) 



Hereafter we call this rate "G-rate" (see Figure [2|). 
The rate is for a cosmological halo forming from 
a very high-cr peak, and is slightl y larger than 



that adopted in a previous work (jSpolvar et al 
20091) . We also run models with three different 
accretion rates 1.0, 0.5, 0.2xG-rate. The inter- 
mediate value , 0.5 x G -rate, is close to that used in 
Spolvar et al. ( 20091) . We use this model to com- 
pare with their result. 

The energy generation rate of DM annihilation 
is given by 



GeV cm"3 s"^ 



(4) 



where lav) is the DM self-annihilation rate in 



units of cm^ 



" ^ , and rriv 



cle mass in units of GeV 



is the DM parti- 
: 10^ electron volt 
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e.g. iBertone et al. ( 2005 )). Notice that effec 



tively only the ratio, {av)/m^, determines the 
net energy generation rate. We vary the DM 
particle mass as a model parameter as = 
1, 10, 20, 50 100 and 200 GeV , while fixing 
(cru) constant. Spolvar et al.l ( 20081) calculate the 
critical transition gas density above which most 
of the DM annihilation energy is absorbed inside 
the core. The critical density is smaller than the 
stellar density with the WIMP mass 1 GeV to 10 
TeV, and thus we assume that the released en- 
ergy from DM annihilation quickly thermalize the 
star. Through the self-annihilation process, DM 
changes into multiple species including neutrinos, 
which escape away from the star and the cloud. 
We assume that about one third of DM annihila- 
tion energy is carried away from the star, and thus 
the effect ive energy generatio n rate is (2/3)(3dm 
(see, e.g. ISpolvar et aD (|2008h ). 

All our calculations start from a protostar with 
M ~ 10 Mq. As the initial condition, we calcu- 
late an equilibrium stellar structure for this mass, 
including DM annihilation energy. In order to 
apply the Blumenthal's method for the adiabatic 
contraction calculation, we need the gas density 
profile outside the star. We use the result in the 
cosmological simulation of lGao et al. I (|2007l) which 
shows the gas density p scales with distance r as 



oc r 



-2.3 



Within an integration time step, the 
accreted gas is added to the stellar mass and DM 
annihilation energy is re-calculated. Because the 
DM density increases extremely within the central 
star due to adiabatic contraction, the DM annihi- 
lation rate at the outside of the star is relatively 
small than inside. Thus we consider the effect of 
DM annihilation inside the star only. We stop 
the calculation when the star contracts to be suf- 
ficiently compact and the stellar surface tempera- 
ture reaches 10^ K. The dark star phase ends well 
before entering this regime. Afterwards the star 
is expected to settle quickly on the main-sequence 
phase, mainly powered by nuclear reactions. In 
this study, we do not include nuclear reactions for 
our dark star models, because the nuclear reac- 
tions are ineffective owing to the low temperature 
inside the star during the dark star phase. 

In the next section, we first show the results 
from the run with a fiducial set of parameter val- 
ues. We call this model "base model". The pa- 
rameters of our base model are a power-law gas 



mass accretion rate l.OxG-rate, = 100 GeV 
and {av) = 3 x 10~^^ cm'^ s~^. We will also show 
the results of models which adopt different values 
for these two parameters. 

3. Results 

3.1. Base model 

Figure [3] shows the evolution of the dark star in 
the Hertzsprung-Russell (HR) diagram. The ac- 
creting star grows constantly in mass. In the HR 
diagram, the star moves upward and then turns 
to the left (higher temperature). The solid line is 
for our dark star model, whereas the dashed line is 
for a standard Population III model (without DM 



annih ilation energy; no-DM model, lUmeda et al 
(|2009 f)). Symbols indicate the times when the stel- 
lar mass is M = 100, 200, 600, 800 and 1000 Mq. 
Clearly the dark star model moves on a very differ- 
ent path from that of the no-DM model, showing 
a significant effect of the DM annihilation energy. 
In the no-DM model, the star contracts gravita- 
tionally and becomes hot quickly until the first 
symbol (100 M0). The dark star, on the other 
hand, has essentially the same effective tempera- 
ture because DM annihilation energy prevents the 
star from gravitational contraction. When the no- 
DM star reaches ^ 100 M©, it expands, and turns 
to the lower temperature side in the HR diagram. 
There, the hydrogen burning becomes effective to 
supply energy to stop the contraction and the star 
expands slightly. Subsequently the star lands on 
the main-sequence phase. 

Contrastingly, the dark star remains "cool" 
while it is powered by DM annihilation. This is 
because the DM annihilation energy (Eq. U) is in- 
dependent of the stellar temperature. Even if the 
stellar interior remains at low temperatures, DM 
annihilation can produce enough energy to sustain 
the stellar structure. As long as this condition is 
met, the star cannot contract while still grows in 
mass by accretion. Interestingly, by the time when 
the dark star reaches main-sequence, the stellar 
mass is already '-^ 1000 Mq, which is much larger 
than the no-DM case (100 - 200 Mq). 

Figure 2] shows the evolution of the DM mass 
inside the star and its time derivative. The left 
panel shows the total DM mass inside the star, 
whereas the right panel shows the ratio of the rate 
of DM annihilation in mass to the rate of DM mass 
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increase (within the star) owing to adiabatic con- 
traction: 



AM 



DM 



dt 



dM, 



DM 



annihilation/ 



dt 



(5) 



contraction 



In Figure m we show the resuhs for an additional 
model in which the depletion of DM in the star is 
not taken into account. Although this case might 
appear unrealistic, it is useful to see when DM 
depletion becomes substantial. The solid line rep- 
resents our "base model" while the dashed line 
shows the result of the no-depletion case. Up to 
M - 200 - 300 Mq, there is almost no differ- 
ence between the two cases. At M > SOOMq, 
however, DM depletion becomes significant and 
then the total DM mass actually starts decreas- 
ing. The evolution thereafter is interesting. The 
DM "fuel" inside the star runs short to sustain 
the star. The star stops expanding, the gas den- 
sity increases more efficiently, and then DM den- 
sity also increases again. Then DM annihilation 
consumes rapidly the DM fuel inside the star. Fi- 
nally, the DM annihilation energy cannot sustain 
the star, marking the end of the dark star phase. 
The star collapses and will eventually reach the 
main-sequence phase. 

Figure [5] shows the radial profiles for various 
quantities of gas {left) and DM {right) when the 
stellar mass is M — 200 (solid lines), 800 (dashed 
lines) and 1000 M© (dotted lines) respectively. 
The horizontal axis is the radial distance from the 
stellar center divided by the stellar radius. 

The top panels show the evolution of gas den- 
sity and DM density. Between M = 200 to 
800 Mq, the star has an extended structure and 
the central density does not increase much. The 
DM density increases by adiabatic contraction but 
actually decreases slightly at the inner most part 
owing to annihilation. At the final contraction 
phase {M = 800 to 1000 M©), both the gas den- 
sity and the DM density increase substantially. 

The middle panels show the evolution of en- 
closed mass of gas and DM. Although the gas mass 
profile stays roughly unchanged, the DM mass de- 
creases dramatically during the final phase from 
M = 800 to 1000 Mq. Note that the horizontal 
axis in the plots shows a normalized radius R/R^. 
The stellar radius i?, itself changes significantly 
over the plotted range of evolutionary stages. Be- 
cause the dark star phase ends with the runaway 



burning of DM inside of the star, the total amount 
of DM when M = 1000 Mq is already very small. 

The bottom-left panel shows the DM annihi- 
lation rate and the bottom-right panel shows the 
total luminosity generated by the DM annihilation 
within the radius. These panels show the energy 
generating efficiency from the DM annihilation. 
Again, we see that the annihilation rate is low in- 
side the star at the final stage M — 1000 Mq and 
the enclosed DM luminosity at the stellar surface 
(at logioRadius/R* — 0) is smaller than in the ear- 
lier phases. The stellar mass increases by gas ac- 
cretion but the DM energy supply decreases; this 
causes the star to contract. 

Figure [5] shows the evolution of some basic stel- 
lar quantities which characterize the dark star. We 
compare the results of our base model (solid lines) 
with those of no-DM model (dashed lines) and of 
no-depletion model (dotted lines). Until the star 
grows to M ~ 200 — 300 M© , the base model and 
the no-depletion model appear very similar, show- 
ing again that DM depletion is negligible in early 
phases. After the star grows to ~ SOOMq, we see 
small but appreciable differences between the base 
model and the no-depletion model. As has been 
shown in Figure |4l DM is consumed by annihila- 
tion inside the star whereas the supply by adia- 
batic contraction is slow owing to the small gas 
mass accretion rate. As the star becomes more 
massive, it needs more DM to produce necessary 
energy to sustain gravitational equilibrium. When 
the DM supply becomes insufficient, this quasi- 
stable dark star phase cannot be sustained. It oc- 
curs at M ~ 600 Mq for the base model. The 
stellar properties, however, do not change imme- 
diately until up to M ~ 900 Mq. At around 
M > 900 Mq, the DM inside the star burns out 
rapidly, the star begins to collapse, and the central 
temperature increases rapidly. Finally the central 
temperature reaches 10^ — 10® K, the nuclear burn- 
ing will soon start, and the star will eventually 
land on the main-sequence. 

Table [T] summarizes the basic stellar proper- 
ties at several characteristic phases. Note that the 
elapsed time difference (the right column) is large 
between 800 and 1000 Mq. The dark star phase 
continues for about ~ 0.2 Myr. Throughout the 
run, the DM mass is very small compared with the 
stellar gas mass. Dark stars can be sustained by 
the DM of only 0.1% of their total mass Af,. 
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3.2. The effect of gas mass accretion rate 

We now explore a few parameter spaces. Our 
dark star model has essentially two parameters, 
the gas mass accretion rate and the DM particle 
mass. First, we examine the effect of the gas mass 
accretion rate. We adopt three accretion rates, 
1.0, 0.5, and 0.2xG-rate (Figure d]). "G-rate" 
is the gas accretion rate in our base model (Eq. 
[3]). In general, less DM is attracted toward the 
star for smaller gas accretion rates. We consider 
only smaller accretion rates than G-rate because 
it seems reasonable to assume that various proto- 
stellar feedback effects, such as ionizing radiation 
from the star, can only reduce the gas mass accre- 
tion rate. 

Figure [7] shows the stellar evolution for the 
three models with different accretion rates in the 
HR diagram. Initially, the stars are puffy, cool and 
evolve along virtually the same path. Tables [5] and 
[3] show the basic stellar properties when the mod- 
els with the reduced accretion rates reach the end 
of the stable phase where T^s = 10^ K, and where 
Tcs = 10^ K, respectively. The latter phase is the 
end of our calculations. As is expected naively, the 
final stellar mass is smaller for lower gas accretion 
rates, but the period of the dark star phase is actu- 
ally longer. In the reduced accretion models, the 
star needs a longer time to reach the same stel- 
lar mass. Then the DM density inside the star is 
smaller because the gas contraction is slow. Con- 
sequently, the star with a smaller gas accretion 
rate ends the dark star phase and begins to con- 
tract when its mass is smaller than in the base 
model. 

Figures [5] and IH] show clearly the effect of re- 
ducing the accretion rate. There are essentially 
no differences among the three models when the 
stellar mass is small. The DM fuel inside the star 
increases very similarly (Fig. [8^) because the rate 
of DM decrease to DM increase is very small (Fig. 
[8)3) . For the smallest accretion rate, the dark star 
ends when the stellar mass reaches ^ 500 Mq. 
Overall, the lower the accretion rate is, the smaller 
the DM mass is inside the star. The dark star 
phase lasts longer but ends at lower masses. 

3.3. The effect of DM particle mass 

Finally, we study models with different DM par- 
ticle masses. We adopt the base accretion rate of 



1. Ox G-rate which is larger than those adopted in 
Spolvar et al.l ( 20091 ) (see Figure [2|). so our results 



are slightly different from theirs, especially for a 
small DM particle mass case. We discuss this issue 
later in this section. 

We run six models with = 1, 10, 20, 50, 100 
and 200 GeV. Note that, from Eq. H the DM 
annihilation energy is inversely proportional to 
. Thus the characteristic features of dark stars 
appear strongly in runs with small DM masses. 
Figures [lOl [TT] and [12] show the results for the 
five models. For smaller DM particle masses, the 
DM annihilation rate is large, and thus the star 
expands more and the temperature and density in 
the star are lower. The period of the dark star 
phase is also prolonged for smaller m^, and the 
final stellar mass gets larger than our base model. 

Tables |4] and [5] list the basic properties of the 
stars when the surface temperature Tcff reaches 
lO'' K and 10^ K, respectively. The final stellar 
mass differs by about a factor of two among 
the five models; = 776 to 1370 GeV in Ta- 
ble m The overall trend is as describe above, 
but the lightest particle model is worth describ- 
ing more in detail. When the star enters main- 
sequence, the stellar mass in = 1 GeV model 
is similar to that of 10 GeV model, although it 
has 10 times larger DM energy generation rate. 
The reason is that the DM burning rate is too 
large in = 1 GeV model. When the DM con- 
sumption rate exceeds the DM supply by adia- 
batic contraction (see discussion in Sec. 3.1), the 
star runs out of the DM fuel and begins to con- 
tract. The DM consumption becomes even more 
efficient then. This final "runaway" phase takes 
place faster in 1 GeV model than in 10 GeV model. 
Consequently, the star lands on the main-sequence 
early in 1 GeV model (see the last column of Table 
[5] which gives the elapsed time). 

Interestingly, the — 1 GeV model predicts 
the most luminous star, reaching L ^ some 10^ Lq 
at the peak (Fig. [T2]) . The extremely high lumi- 
nosity is, however, still smaller than the upper- 
limit given by the Eddington Luminosity Ledd^ 



AttcGAL 



^EDD 



(6) 



where the opacity k determines essentially the crit- 
ical luminosity. Figure [T3] shows the ratio of the 
energy generation rate of DM annihilation ("lu- 
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minosity") to the Eddington luminosity for our 
five models with different DM particle masses. 
The ratio stays always less than unity, even for 
1 GeV model. During the dark star phase (Toff < 
lOjOOOK), the opacity k at lower temperatures is 
mainly contributed hy H~ ion which has a small 
value. 

4. Discussion 

4.1. Dark Star Models 

We compare our results with tho se of previous 



works . Unlike previous studies fe.g.. lSpolvar et al 
we solve full radiative transfer and follow 



self-consistently the stellar structure and its evo- 
lution. A direct comp arison can be made w ith the 
results of Table 2 in 'Sp olvar et all (|2009t ) (here- 
after "SP09"), by using the result of our 0.5 xG- 
rate model (see Sec. 3.2; hereafter "05G"). The 
other model parameters are set to be the same, 
= 100 GeV and {av) = 3 x IQ-^^ cm^ s-\ ex- 
cept for the gas mass accretion rate. The gas mass 
accretion rates in these two models (Figure [2]) are 
almost the same at M^, > 100 Mq, which is the 
main period of the dark star phase. The results 
of 05G and SP09 are roughly similar; in stellar 
radius, luminosity and other quantities, the differ- 
ences are less than 20% at Af* = 106 Mq (the first 
line in SP09), and remain within a factor of two 
at later stages. During calculation, 05G model 
has a slightly expanded structure and is some- 
what cooler than SP09 model. On the other hand, 
the difference of the central temperature is rela- 
tively large between the two models. The value 
in 05G becomes lower than in SP09 especially at 
the final stages; at = 716, 756, 779 Mq, 
Tccn = 2.4, 9.3, 132 (106)K in 05G whereas 
Teen = 15, 78, 280 (lO^)K in SP09. Because of 
the lower central temperature, nuclear burning is 
unimportant in our calculations. 

4.2. Possible Range of WIMP DM mass 

The WIMP mass significantly affect the dura- 
tion of the dark star phase because the DM an- 
nihilation energy is inversely proportional to the 
WIMP mass m^, as can be seen in Eq. S) We 
have already shown the effect in Section 3.3 for 
a typical DM mass range, from 1 to 200 GeV. 
However, in supersymmetric particle physics mod- 
els, the WIMP DM mass may be 10 TeV or even 



greater. To estimate the evolutionary path of dark 
star models with heavier WIMPs, we have run our 
no-depletion model for two larger WIMP masses, 
= 1 and 10 TeV. Figure [TIk and [Tib show 
the result. At the beginning of the calculation 
with M ~ 15Mq, all the models are in the dark 
star phase and on the Hayashi-line in the H-R di- 
agram. The left panel of Fig. [M] shows that, at 
this phase, all models have Teff — lO'^ '' K. For 
= 10 TeV, the Pop. Ill star leaves the dark star 
phase early. Because annihilation rate is small for 
a large WIMP mass, the stellar mass at the end of 
the dark star phase remains small. 

All the dark star models are on Toff — lO'^'^ K 
line at first (see also Fig. [TU|) and leave the 
line gradually as mass- accretion continues. To 
quantify the evolution, we define two charac- 
teristic stellar masses, M3.8 and M4.0, the stel- 
lar masses at which the surface temperature 
reaches T^s ~ lO^-^ and 10^ ° K. The former 
indicate the phase where the star deviates from 
the Hayashi-line. The later indicates when a 
star ends the stable dark star phase. The re- 
sult is, M3.8 = (1145, 821, 402, 74, 26) Mq 
and M4.0 = (1349, 1325, 821, 270, 50) Mq for 
= (1> 10, 100 GeV, 1, 10 TeV), respectively. 
We conclude that Pop. HI stars with a large WIMP 
mass up to 10 TeV can path through the dark star 
phase, although the duration of this phase are very 
short. The final stellar masses are small for the 
dark stars with most massive WIMPs. 

Note that we have used no-depletion models 
here for simplicity, because dark matter depletion 
effect is small for the early stages of the dark star 
evolution. If we take the DM depletion into ac- 
count, a dark star ends the stable phase earlier 
and above two quantities A/3.8 and M^^q become 
smaller. 

5. Conclusion 

We have studied the pre-main-sequence evolu- 
tion of dark stars by following the self-consistent 
stellar evolution. We have suitably modified the 
code to incorporate the energy generation from 
spherically distributed DM. 

Our base model with = 100 GeV and 
dM/dT = l.Ox G-rate Mq yrs^^ shows the char- 
acteristic features of the dark star phase; the large 
DM annihilation energy expands the star, making 
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the star to be in gravitational equilibrium. The 
lower temperature is one of the peculiar proper- 
ties of the dark star. This stable phase continues 
until the energy supply from DM annihilation be- 
comes insufficient to maintain the stable structure. 
Finally, the star collapses rapidly and reaches the 
main-sequence phase. At this point, the stellar 
mass has grown up to M - 900 - 1000 Mq. 
Such features in dark star phase are a ll consistent 
with t h e findings of p r eviou s works of llocco et al 



(|2008l ): ISpolvar et al] (|2009( ) who employed much 
simpler stellar models. 

The dark star model has effectively two param- 
eters. One of them is the gas mass accretion rate 
which determines the evolution of the gas and 
DM distribution. Cosmological simulations pre- 
dict a variety of gas accretion rates. For a small 
accretion rate, the period of the dark star phase 
increases whereas the final stellar mass decreases 
M ^ 500 Mq. Note however that the mass is still 
larger than the standard Pop III (no-DM) case 
(100 - 200 Mq). Another parameter is the DM 
particle mass which determines the energy gen- 
eration rate. For a small DM particle mass, the 
period of dark star phase becomes longer and the 
final stellar mass becomes larger (see Table [S]). 

All models pass through the dark star phase 
and this phase is maintained by a little DM fuel 
which is less than 0.1 % of the stellar mass. If 
the first stars in the universe have undergone such 
phase, there are exotic stars which are cool and 
massive in the early universe. 

Formation of very massive dark stars has an 
important implication. Such stars eventually col- 
lapse gravitationally, to form massive black holes 
with masses as large as 1000 Mq. The rem- 
nant black holes can also grow by accretion or 
by mergers to seed super-massive black holes. 
The exotic feature of dark star's appearance, lu- 
minous and cool, may be able to use as the 
powerful clue to search such stars. Our cal- 
culations show that first stars can grow to be 
dark stars with luminosity ~ few 10^ Lq at 
most and they will no t be d etectable by JWST. 
However, iFreese et al. 1 (|20in[) estimate for super- 
massive dark stars (SMDSs) and such very mas- 
sive and bright st ars (about 10^~^^ £q ) can be de- 
tected by JWST. ISandick et all (|2010[ ) argue that 
dark star remnants might survive to the present 
day in the Milky Way, leaving 7-ray signatures 



from DM annihilation. 

In future work, we will include the nuclear re- 
action to calculate the first star evolution from 
the dark star phase to the main-sequence phase 
completely. It would be also interesting to include 
capture of DM particles. Once the stars have con- 
tracted and stop growing in mass, there will be 
no more DM supplied by adiabatic contraction. 
However, DM particle may still be captured by 
the star. Previous studies showed that this pro- 
cess can make the star back in the dark star phase 
again. We will explore the effect of DM capture 
and make a complete evolutionary model of the 
first stars with DM annihilation. 
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the MEXT of Japan and the Grants-in-Aid for 
Young Scientists (S) 20674003 by the Japan So- 
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from World Premier International Research Cen- 
ter Initiative (WPI Initiative), MEXT, Japan. 
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Fig. 1. — Dark matter density profile calculated by an adiabatic contraction model. The solid and dashed 
lines are calculated directly from our cosmological iV-bod y simulation, whereas the dotted line is calculated 
from the initial profile data using the analytical method of lBlumenthal et"ar (1986). The dashed and dotted 
profiles agree remarkably well. 
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Fig. 2. — Gas mass accretion rates adopted in our stellar evolution calculations. "G-rate" is the accretion 
rate of Eg. El We ad opt three rates 1.0, 0.5 and 0.2x G-rate. The other two lines are the rates adopted in 
Spolvar et al.l ( 20091) . for comparison. 
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Fig. 3. — Hertzsprung-Russell (HR) diagram for our base model. Symbols indicate the stellar mass of 
M = 100, 200, 600, 800 and 1000 Mq. The sohd hne is for the d ark star mode l ("bas e model"), whereas 
the dashed line is for the standard Pop III model (no-DM model of lUmeda et al. I (t2009h ). 



Base model 




Mass [Mq] Mass [MJ 

Fig. 4. — DM mass inside the star (left panel) and the DM consumption rate by the annihilation (right 
panel) as a function of stellar mass. We normalize the DM consumption rate dM z) m /dt\ ^annihilation by the 
DM mass supply by adiabatic contraction dA/DM/c^^l contraction- The solid line is for the dark star model 
("base model"), whereas the dotted line is for no-depletion model. 
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Fig. 5. — Radial profiles of various quantities for our base model when the stellar masses are M = 200, 800 
and 1000 M©. In all the panels, the horizontal axis shows the ratio of radial distance normalized to the 
stellar radius (see Table [T]). Top panels show the gas and DM density profiles. The middle panels show 
the gas and DM enclosed mass. The bottom panels show DM annihilation rate and DM luminosity from 
annihilation. 



13 



Base model 

No-DM model 

no-depletion 





4 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




3 






2 






1 






3 


l"m|mm|mm|mm|mm|m 




1 




cm 


-1 




§ 


-3 




Q. 


-5 






-7 





200 400 600 8001000 
Mass [M ] 



' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' 




il 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




200 400 600 8001000 
Mass [M ] 



Fig. 6. — Evolution of basic stellar quantities. The horizontal axis indicates the stellar mass. The vertical 
axis shows, respectively, the stellar radius, the central gas density, the luminosity , the surface temperature 
and the central temperature in the logarithm scales. The solid line is for the base model, the dashed line is 
for no-DM model and the dotted line is for no-depletion model. 
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Fig. 7. — HR diagram for models with three accretion rates, 1.0, 0.5 and 0.2x G-rate, shown in Figure[2] 
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Fig. 8. — DM mass and the DM consumption rate by the annihilation as a function of stellar mass for three 
accretion rates models, as in Figure ID 
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Fig. 9. — Evolution of basic stellar quantities for dark star models with reduced accretion rates. AH vertical 
axes are plotted in logarithm scale. The plotted quantities are the same as in Figured 
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Fig. 10. — HR diagram for models with variation of DM particle masses. Each lines show models with 
DM particle masses — 1, 10, 20, 50, 100, 200 GcV and the no-DM model (the same as in Figure [3]), 
respectively. 




Fig. 11. — DM mass and the DM consumption rate by the annihilation for the variation of DM particle 
masses, as the same plot of Figure ID 
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Fig. 12. — Evolution of stellar values for dark star models with parameters of DM particle masses. All 
vertical axes are plotted in logarithm scale. As for Figure [51 but for different DM particle masses. 
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Fig. 13. — We plot the ratio of DM annihilation luminosity to Eddington luminosity (left) and the opacity 
at photosphere as a function of stellar mass (right) for different DM particle masses. In all the models, the 
luminosity ratio stays less than unity. 
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Fig. 14. — Stellar values of no-depletion models with large WIMP mass = 100 GeV, 1 and 10 TeV. 
Left panel shows the evolutions of the stellar surface temperature and right panel shows ones of the stellar 
radius. All models are on the dark star phases at first. 
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Table 1 

Stellar properties for the "base model" calculation 



[Mq] 


R* [R-©] 


L* [Lq] 


Tcff [K] 


Teen [K] 


Mdm [Mq] 


t [Myr] 


200 


1.51E3 


1.83E6 


5.49E3 


4.69E5 


0.166 


0.017 


600 


1.51E3 


6.98E6 


7.67E3 


9.33E5 


0.313 


0.097 


800 


1.19E3 


1.08E7 


9.65E3 


1.42E6 


0.287 


0.153 


1000 


21.98 


3.24E7 


9.33E4 


9.25E7 


3.75E-4 


0.218 



Note. — We tabulate the main results for the base dark star model at stellar 
masses, M = 200, 600, 800 and 1000 M©. The columns show, from left to right, 
the stellar radius, luminosity, surface temperature, central temperature, DM mass 
inside the star, and the time elapsed from the beginning of the calculation. 



Table 2 

Stellar properties for three accretion models at T^s = 10^ K 



Accretion Rate 


[Mq] 


R* [R©] 


[Lq] Tefr [K] 


Teen [K] 


Mum [Mq] 


t [Myr] 


1.0 X G-rate 


821 


1.14E3 


1.15E7 1.00E4 


1.51E6 


0.277 


0.160 


0.5 X G-rate 


682 


9.87E2 


8.67E6 1.00E4 


1.56E6 


0.205 


0.237 


0.2 X G-rate 


492 


7.87E2 


5.55E6 1.00E4 


1.63E6 


0.125 


0.351 



Note. — Stellar properties for three accretion models when the stellar surface temperature reaches 
Teff = 10^ K. After this phase, the star begins to run out the DM fuel and begins to gravitationally 
contract. 



Table 3 

Stellar properties for three accretion models at Teff = 10° K 



Accretion Rate 


[Mq] 


R* [R-o] 


[Lq] Teff [K] 


Teen [K] 


MdM [Mq] 


t [Myr] 


1.0 X G-rate 


1002 


19.27 


3.32E7 1.00E5 


1.07E8 


2.43E-4 


0.220 


0.5 X G-rate 


775 


16.55 


2.43E7 1.00E5 


1.08E8 


1.46E-4 


0.291 


0.2 X G-rate 


534 


13.44 


1.62E7 1.00E5 


1.09E8 


7.65E-5 


0.401 



Note. Stellar properties for three accretion models when the stellar surface temperature reaches 
Teff = 10^ K. By this phase, the star has contracted sufhciently and the density and temperature 
has risen up enough to start the hydrogen burning. In this point, the star is not supported by DM 
annihilation energy, so we stop the calculation. 
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Table 4 

Stellar properties for some DM mass models at T^s = 10^ K 



[GeV] 


[Mq] 




[Lq] 


Teff [K] 


Teen [K] 


Mdm [Mq] 


t [Myr] 


1 


1349 


1.82E3 


2.92E7 


1.00E4 


1.11E7 


0.096 


0.354 


10 


1325 


1.59E3 


2.23E7 


1.00E4 


1.33E6 


0.254 


0.344 


20 


1214 


1.51E3 


2.02E7 


1.00E4 


1.42E6 


0.297 


0.298 


50 


1009 


1.32E3 


1.55E7 


1.00E4 


1.46E6 


0.307 


0.222 


100 


821 


1.14E3 


1.15E7 


1.00E4 


1.51E6 


0.277 


0.160 


200 


536 


8.76E2 


6.81E6 


1.00E4 


1.47E6 


0.198 


0.081 



Note. — Stellar properties for DM mass variation models when stellar surface temperature 
reaches Tes = 10* K; the end of the stable dark star phase and the start of transformation into 
main-sequence star. 



Table 5 

Stellar properties for some DM mass models at T,-.[[ = 10"'' K 



[GeV] 


[Mq] 


-R* [R©] 


L* [Lq] 


Teff [K] 


Teen [K] 


Mdm [Mq] 


t [Myr] 


1 


1370 


23.86 


5.03E7 


1.00E5 


2.49E8 


3.47E-6 


0.362 


10 


1381 


23.32 


4.83E7 


1.00E5 


1.07E8 


4.96E-5 


0.367 


20 


1295 


22.45 


4.48E7 


1.00E5 


1.07E8 


8.92E-5 


0.333 


50 


1138 


20.85 


3.86E7 


1.00E5 


1.07E8 


1.59E-4 


0.269 


100 


1002 


19.27 


3.32E7 


1.00E5 


1.07E8 


2.43E-4 


0.220 


200 


776 


18.91 


3.16E7 


1.00E5 


9.95E7 


3.88E-4 


0.146 



Note. — Stellar properties for DM mass variation models when stellar surface temperature 
reaches T^g = 10^ K; the end of calculation. 
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